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We report on the spontaneous production of fluxons in annular Josephson tunnel junctions during 
a thermal quench in the presence of a symmetry-breaking magnetic field. The dependence on field 
intensity B of the probability /i to trap a single defect during the N-S phase transition depends 
drastically on the sample circumferences. We show that this can be understood in the framework 
of the same picture of spontaneous defect formation that leads to the experimentally well attested 
scaling behaviour of /i with quench rate in the absence of an external field. 

PACS numbers: 03.70.+k, 05.70.Fh, 03.65.Yz 



I. SPONTANEOUS FLUXON PRODUCTION IN 
AN EXTERNAL FIELD 

Some time ago it was proposed by Kibble fH and Zurek 
J?] that the domain structure after a continuous phase 
transition is determined by its causal horizons. Since 
then there has been a sequence of experiments in con- 
densed matter systems to attempt to confirm this by 
counting the topological defects produced at a quench, 
whose number can be correlated to phase boundaries. 
Such confirmation lies in the predicted scaling behaviour 
of the average defect separation ^ at the time of their pro- 
duction as a function of the quench time (inverse quench 
rate at the critical temperature Tc) tq. Supposing that 
the 'equilibrium' correlation length £,eq{T) of the order- 
parameter field, and its relaxation time t{T), diverge at 
T = Tr as: 
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where the scaling exponent cr > belongs to universality 
classes determined by the (adiabatic) critical exponents 
u and 7, and v, 7, ^0 and tq depend on the detailed 
microscopic behaviour of the system. The coefficient of 
proportionality in ([1]) is an efficiency factor which varies 
with the system, from a few percent for high temperature 
superconductors [1,0] to full efficiency for superfluid '^He 

In the last sever al y ears we have performed a set of 
experiments [1, H, [13, [HI on planar annular Joseph- 
son Tunnel Junctions (JTJs) to test the scaling law ([1]). 
Specifically, a planar Josephson junction comprises two 
superconducting films separated by an insulating oxide 
layer. We assume continuity in the density of Cooper 
pairs across the oxide, but allow for a discontinuity (j) 
in the phase of the effective order parameter field. Once 
the transition is completed the lossless Josephson current 
density is J = Jc sin (j), for critical current density Jc. For 
such JTJs the defects are fluxons (or antifluxons) corre- 
sponding to a change A(j} = ±2m: in along the oxide 
layer |12l |. The integer ±n is the so-called winding num- 
ber. The Swihart velocity provides the requisite causal 
horizons [l^ . 

Most simply, for small annuli of circumference C <C f , 
the trapping probability /i for finding a fluxon (or /_i 
of finding an antifluxon) is taken to be: 



/i = ./-! = C/e'cx (C/6)(TQ/ro)- 



(2) 



In our experiments we measure /i — /i -I- the like- 
lihood of seeing one fluxon or antifluxon. We have car- 
ried out statistical measurements to determine the de- 
pendence of /i on TQ for a large number of high qual- 
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FIG. 1: Log-log plot of the frequency /i of trapping single 
fiuxons versus the quenching time tq for an AJTJ of circum- 
ference 0.5mm (closed circles) and for an AJTJ with circum- 
ference 1.5mm. The best fits through the data (solid and 
dashed lines) show that scaling with a = 0.5 is totally robust. 
Both samples had equal critical current density Jc- 



FIG. 2: Dependence of the single fluxon trapping frequency /i 
for quench time tq = 5s for an AJTJ of circumference C — 
0.5mm in the presence of a magnetic field B perpendicular 
to the barrier plane. The solid line corresponds to the case 
iV = 1 in fi{N,n) of Fig.5, as derived from Eg. ((22)) . 



ity AJTJs with circumferences varying from 0.5 mm to 
3.14: mm. All samples had equal critical current density 
Jc{T = 0) ~ 60 A/ cm^ (corresponding to a Josephson 
penetration depth Xj{T = 0) ~ SO/im), yielding the same 
^0 and To . The quenching time tq could be changed over 
several orders of magnitudes from tens of milliseconds to 
tens of seconds, using the methods described in [lol. 

In Fig. 1 we show /i(tq) for annuli of radius 0.5mm 
(lower plot) and 1.5mm (upper plot). The results for the 
larger annulus have not been shown before. There is no 
doubt that, for small C/^, the probability /i of finding 
a single fluxon shows scaling behaviour of the type 
with a = 0.5 to high accuracy, some of whose details 
were given m [13, OM- This value of a is as we would 
expect for realistic junctions for which the fabrication 
leads to a proximity effect [3, [HI and for which the Swi- 
hart velocity does not show complete slowing down. The 
efficiency factor is approximately unity for the smaller 
annulus. More details are given in [ll| . 
What concerns us in this paper is how fiuxons form in the 
presence of an externally applied magnetic field B that 
explicitly breaks the symmetry of the theory, whereby 
/i 7^ This is a crucial ingredient in the analy- 

sis of unbiased fluxon production in JTJs because, de- 
spite our best efforts, we cannot preclude the possibility 
of stray magnetic fields in the experimental equipment 
(e.g. a magnetised screw or an incomplete shielding of 
the earth's magnetic field). In fact, in presenting the data 
in Fig. 1 we have taken the effects of static stray fields em- 
pirically into account, by applying an external field until 
such stray fields are neutralised. It is only then that we 
obtain In all cases the external magnetic field B 

was applied perpendicular to the junction plane. This 



choice of field orientation is mainly due to the fact that a 
transverse field (due to demagnetization effects) is more 
effective, by almost two orders of magnitude, than an in- 
plane field in modulating the junction critical current Ic 
jig and trapping frequency. Furthermore, under partic- 
ular conditions [Tg a transverse magnetic field allows to 
discriminate between fiuxons and antifluxons. 

As a result, we have built up a substantial collection of 
data showing the dependence of /i on both t and B, that 
we shall discuss in the remaining sections. In particular, 
we shall concentrate on two representative datasets that 
refer to high quality Nb/Alox/Nb - Nb AJTJs quenched 
at the same quench rate (tq — 5s), but having different 
circumferences, i.e., C = 0.5mm and C = 2.0mm, shown 
in Fig. 2 and Fig. 3, respectively. Details of the samples' 
electrical and geometrical parameters and of the experi- 
mental setup can be found in |11|. 

We observe that the increase in circumference has a 
dramatic effect on the single trapping frequency /i. For 
the 0.5mm long AJTJ we find a central minimum with 
two side-peaks. Such a double-peaked dataset has been 
used for any single data point in Fig.l, with fi being 
read off from the central minimum, typically displaced 
slightly from B = because of the aforementioned stray 
fields in the equipment (typically several tens of nT). 

On the other hand, for the 2.0mm long sample (not 
represented in Fig.l) we only have a central peak. Super- 
ficially this is strange, since it shows that the probability 
of seeing a single fluxon decreases as the external field 
increases, but we shall understand this as a consequence 
of an increased ability to create more than one fluxon. 
We note that the values of the magnetic field required to 
change /i significantly are very small when compared to 
the field values needed to modulate the Josephson cur- 
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FIG. 3: Dependence of /i for quench time tq = 5s for an 
AJTJ of circumference C = 2.0mm in the presence of a mag- 
netic field B perpendicular to the barrier plane. The solid 
line corresponds to the case TV = 16 in /i(7V, n) of Fig.5, as 
derived from Eq. (f22)l . 

rent Ic of the samples, whose first minimum occurs at 
field values of several fiT. Further, the larger the ring 
size the larger is the efi'ect of a given magnetic field. 

In the remainder of this paper we shall show how the 
results of Figs. 2 and 3 can be understood, and qualita- 
tively predicted, in a framework that implies the scaling 
behavior of Fig.l. As such, these results are mutually 
supportive but require an alternative formulation of the 
KZ scenario, which we now provide. 

II. CAUSALITY VS. INSTABILITY 

At first sight the results for fluxon production in an 
external field have little or nothing to do with the Kibble- 
Zurek scenario. However, the original scenario of bound- 
ing domains by causal horizons is just one way of saying 
that, qualitatively, systems change as fast as is possi- 
ble. From a different viewpoint we know that continuous 
transitions (like the one here) proceed by the exponen- 
tial growth of the amplitudes of unstable long- wavelength 
modes of the system. This growth is strongly suppressed 
by self-interaction once the system is close to the ground 
states of its symmetry-broken phase. 

Exponential growth is as fast as it gets, corresponding 
to linearising the equations of motion. Thus, provided 
there is enough time for such rapid growth before back- 
reaction (self-interaction) stops it, we can understand 
how systems can change as fast as is possible without 
invoking causal horizons directly. There is a corollary to 
this. The linear behaviour of the order parameter fields, 
while the transition is taking place, follows from their 
behaving as Gaussian random variables. In fact, for an 
idealised situation in which amplitude growth is stopped 



by implementing a rapid back-reaction that chokes off 
growth instantaneously, it can be shown [3 [S that 
the assumption of Gaussianity leads to scaling behaviour 
like that of ([1]), with the same values of a as would be 
obtained from causal bounds. What may seem surprising 
is that, even when back-reaction is included more realis- 
tically, numerical simulations [2l|, [2^ show that the 
behaviour is still essentially the same. The a exponents 
of causal reasoning are recovered. There are, however, 
two major differences between fastest amplitude growth 
and causal bounds. The growth starts after the transition 
has begun, whereas causal bounds can be imposed both 
before or after the transition has begun, usually to the 
same effect 0, Q . Numerical simulations show [l^, HI] , 
without a doubt, that it is the behaviour of the system af- 
ter the transition has begun that determines the domain 
structure. This is necessary in our context since, with no 
Josephson effect for T > Tc, we could not have invoked 
causality before the transition. Further, the assumption 
of Gaussianity gives us more than causality, in that it 
reintroduces the role of the Ginsburg temperature, at 
which thermal fluctuations become important, where ap- 
propriate. This provides one natural explanation for the 
failure to observe vortices in quenches of "^He |19|] while, 
by a similar argument, permitting spontaneous vortex 
production in ^He [1, Q and superconductors 

However, since equations of motion are, by construc- 
tion, causal, these viewpoints are largely complementary 
once these caveats are taken into account. When appro- 
priate we will invoke both mechanisms in our subsequent 
discussion. 

The simulations cited above are largely for systems 
with global symmetry breaking, simpler than supercon- 
ductors. In fact, although local breaking gives a very dif- 
ferent domain structure, the idea of the transition being 
driven by instabilities survives. Again, the exponents are 
those of causal arguments [l^l except that there is a ad- 
ditional mechanism[25| for the spontaneous production 
of flux in which magnetic fleld just freezes in by itself, 
with behaviour very different to that of ([T]). For annu- 
lar JTJs this further mechanism does not arise because 
of the thinness of the oxide layer through which fluxons 
protrude and the scaling behaviour of ([1]) and (HI) is a 
clean prediction. 



III. BEYOND THE LINEAR REGIME {B = 0) 

Suppose that there is no external symmetry-breaking 
fleld. As long as /i is significantly smaller than unity we 
expect the linear log-plot in tq, as seen in Fig.l. How- 
ever, with /i bounded by unity the linear behaviour will 
soon break down for increasingly fast quenches. Once 
^ — 0(C) the trapping probability of finding net 
fiuxon number m (fiuxons minus antifluxons) increases 
for m > 1, forcing fi to decrease. For the remainder of 
this section we shall extend ^ to predictions for /i and 
higher fm across the whole range of C 
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Since the KZ scenario is appropriate for B = 0, this 
is an ideal testing ground for our two approaches. We 
shall elaborate on these in turn and see that, for many 
purposes, they give almost indistinguishable results. 



A. Independent domains and Gaussian 
probabilities 

In the spirit of the KZ scenario a simple, but informa- 
tive first guess as to how /i and other behave across 
the whole range of C/^ is to divide the annulus into 
N > 2 independent (causal) domains in each of which 
the Josephson phase is a constant. We assume that 
there is no correlation between the values of (f) in adja- 
cent domains but, in calculating the total phase change 
A(j) around the annulus, the geodesic rule is adopted [2a |. 
This means that, when jumping from one domain to the 
other, the shortest path in phase will be taken. The re- 
sult of a quench is then modelled as having the system 
divided up into N domains, each with a randomly chosen 
phase. There is nothing in this ansatz peculiar to JTJs, 
and it is equally applicable to superconductors. As such 
it is an idealisation of a superconducting loop, made out 
of Josephson junctions in series, that has been the object 
of spontaneous flux generation 

Let Ga/(A0) be the probability that the change in 
phase <j> is after M domain boundaries. If the sys- 
tem is made of only two domains, then the lack of phase 
correlation requires: 

Gi(A?!)) -!- for -Tr<A(l)<Tr, 
Ztt 

= for |A0| > TT. 

In fact, the approach of using the KZ picture to set up 
discrete domains in which the order parameter field can 
take random values is one that has been used repeatedly 
for counting defects, at least since its introduction by 
Vachaspati and Vilenkin [27| for counting cosmic strings 
(vortices) in the early universe. This adopts an earlier use 
of causal horizons by Kibble from which [l[ evolved. 
On increasing N, GAr(A(/)) is determined by A'' — 1 self- 
convolutions of Gi: 

Gj\j — Gi * ... * Gi . 

N—l times 

Applying the geodesic rule for the final step in phase 
(from domain A^ back to domain 1, the probability of 
ending with a phase shift of 27rm (i.e. net fluxon number 
m) is: 

d$ Gat ($) = 2^GAr+i(2m7r), (3) 

-7r+2m7r 

where the last equality comes from the definition of Gn- 
To bring this further into correspondence with the KZ 
scenario, we should identify the domain size as compara- 
ble to G/f i.e. G = aN^, where a = 0(1). The value of 



a is not unity, since a discrete domain structure is only a 
crude approximation to a continuous phase at a continu- 
ous transition. Further, the result (P) assumes the causal 
bound is saturated, and we have already commented on 
systems (e.g. high- Tc superconductors) for which the in- 
efficiency of producing flux shows that this is not the case 
0. To take general values of G and f into account we 
need to generalise to non-integer A^. 

GAr(A(/)) already shows rapid convergence to the Gaus- 
sian distribution that arises from the central limit theo- 
rem for A^ > 2. For such A^ the obvious way to proceed 
is to adopt this central limit Gaussian distribution. That 
is, we assume that the total phase change A0 around the 
annulus can be expressed as the sum of a random term 
$ and a geodesic-rule correction (5$. If $ has a normal 
distribution with average $ = and variance ct^(A^) i.e.: 



GAr($) 



^27rcr2(Ar) 



exp - 



/ $2 



(4) 



then a simple calculation enables us to identify (|4]) with 
the central limit distribution of the self-convolutions of 
Gi, as described above, provided that: 

The probability to trap a net number m of defects will 
now be: 



^7r+2m7r 

fm{N) = I d$GAr($). 

-7r+2m7r 



(5) 



Experimentally, in the absence of any external field, 
the most important probabilities are for finding one 
fluxon or one antifluxon. The frequency of no trapping 
fo{N) is: 



foiN) = / Gjv($)d$ = er/( 



), (6) 



and 



f±i{N) = 



37r 



^27rcr2(Ar) 



exp ^ 



$2 



2cr2(iV) 



(7) 



Furthermore, for large a'^ >> say > 20, the trap- 
ping frequencies asymptotically approach zero as: 



f±m{N) 



27r 



(72(A) 



exp - 



27r2m2 



(8) 



Finally, in the same limit, the variance of the discrete 
variable n is: 

aliN) = E,r^„^n2/n(A^) =< >- <7\N)/47r\ (9) 

For A^ < 2 wc require a different approach and turn 
to the consequences of assuming Gaussian stochastic be- 
haviour. 
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B. Gaussian correlations 



If X measures distance along the annulus, (j){x) is peri- 
odic (modulo 27r). The fluxon number density (or wind- 
ing number density) is: 



n{x) = —d^(j){x), 
whereby the net fluxon number n is: 



1 

n= dxn(x) — — Ach, 

Jo ^ ' 27T 



(10) 



(11) 



where A0 is the change in (j). 

For winding number density n {x) the ensemble average 
of the net number of fluxons along an annulus of perime- 
ter C, in the absence of an external field is n = (n) — 0. 

We do not need to adopt any particular form for n (x). 
In the light of our earlier discussion we now assume that 
it is a Gaussian variable until the transition is complete, 
whereby all correlation functions are determined by the 
two-point correlation function {n {x) n [y)) . That is, all 
we shall need for probabilities is: 



{n~) — I dx dy {n (x) n (y)) 
lo 



It follows that: 



\ ' 2P-i(p-l)!^ ' 



(12) 



If fm is the probability of finding net winding number m 
taking both positive and negative values), then: 



(13) 



In order to invert equation p3p , we construct the gener- 
ating function Z(z): 



2\p 



p=0 



[2p)\ 



(n2f)=exp(-z^(n^)/2). 



from (fT2l). On the other hand, from (fTSl 



p—O ^ — oo —oo 

That is, the f„i are the Fourier coefficients of the Gaus- 
sian. In particular, 

h = ^ r dzexp{-z'{n')/2) (14) 



and 



1 f" 

f±i — -z— exp(— z^(n^)/2) C0S2: 

27r 



respectively. 

We observe that, for large (n^), where we can take 
the integration limits to infinity, the trapping probability 
falls off as: 

/±™ « —j=l==exp{-m^/2{n^}). 



In practice, this assumption of a Gaussian stochastic den- 
sity can only be approximate, on two accounts. Less 
significantly, from our earlier comments, it ignores the 
non-linearities of the system. More importantly, it does 
not fully accommodate the periodicity of the annulus, 
to which we shall return later. Nonetheless, it will be 
apparent as to which results are reliable. 

In contrasting Gaussian correlations and Gaussian 
probabilities we see that the definitions are dual to each 
other. In the limit of unrestricted integration (large N) 
they are identical if we identify: 



< a^{N)/ATr^ ^ N/12. 



(15) 



However, for small N there will be differences, as we shall 
see. 

On inserting p5|) . we find that the likelihood /i of 
seeing one fluxon or antifluxon is: 



1 



/i(iV) = - / dz exp(-2^7V/24)cos^ 



(16) 



whereas, assuming Gaussian probability, from ^ we 
find: 



/-(7V) = er/(^)-er/(^). 



(17) 



We know that (fT6|) and (fT7| agree for large N but, as can 
be seen from Figures 4, the qualitative and quantitative 
agreement is striking over the whole range. 

In Fig. 4 we compare /i as a function of N, as given 
by Eq. and Eq. ([T7]) (dashed and dotted lines) respec- 
tively; the dots are the values of /i according to the inde- 
pendent sector model for some integer values of N from 
Eq.Q. Agreement is already good at = 2 and very 
good at TV = 4. 

However, the plots of Fig. 4 are somewhat deceptive for 
small N. This should not worry us since, although the 
analytic expression ^T7\ : 

MN) = 0(VA^)exp(-3/27V), (18) 

vanishes faster than any power, p7|) breaks down there 
by definition. On the contrary, fi{N) of p6|l is linear in 
N for small A'', as we supposed in 

Finally, in the inset of Fig. 4, we also display /o(A^), 
the probability of seeing no flux, for the case of Gaus- 
sian probabilities given in ([6]) . Although we do not show 
it, at the scale of the plot, it is essentially indistinguish- 
able from the result of assuming Gaussian correlations as 
given in (fT4|) . 



fo{N) 



_6_ ^ /tt 



6 



0.5 n 




0.0 I 

10 20 m30 40 50 



FIG. 4: We display the probabilities /i = /i + /_i as a func- 
tion of A*', as given by Ea. (|16|l and Ea. (|17|l (dashed and dotted 
lines) respectively; the dots are the values of /i according to 
the independent sector model for some integer values of A'^ 
from Eq.@. The maximum of /i occurs at N = Nc ~ 10.5, 
for which value /i « 49%. The inset shows the probability 
fo as a function of A*', as given by Eq.(|6]), essentially indistin- 
guishable from that by Ea. (|19p at this scale. 



construction the fm{N) of (fT7|) are automatically positive 
as they must be. 

In a qualitative sense it is of little consequence since, 
throughout the linear regime, the probability /> = 1 — 
/o ~ /ii of seeing more than one fluxon is very small but, 
as a matter of principle, we should impose periodicity to 
render the probabilities positive. We do not have a reli- 
able model in which we can do this but we can get some 
idea by supposing that the phase (j) moved in a double- 
well potential rather than the periodic cos </> potential of 
the sine-Gordon fluxon. In that case, assuming Gaussian 
correlations, the density of defects is proportional [1^ 
to (/"(0)//(0))i/2, where /(x) = (0(a;)0(O)). It is now 
straightforward to impose periodicity, whereupon we find 
behaviour similar to that of (fT8|) for very small N, in that 
it vanishes faster than any power (3Q] . If that were to be 
equally applicable here it is difficult to determine when 
such behaviour might occur, since there is no sign of such 
a collapse in Fig.l, but if this analogy is correct, it will 
repair the minor problem of small negative probabilities 
while leaving the similarity between the two approaches 
at a quantitative level. In particular, as long as we are 
not looking at the small-N behaviour in too much detail, 
as we shall not hereafter, it becomes sensible to use the 
simpler Gaussian probabilities over the whole range. 



We can comfortably use either. 



C. Consequences 

In summary, the assumptions of Gaussian probabili- 
ties (as follows from the KZ picture) and Gaussian cor- 
relations are complementary, with the latter providing a 
(linear) interpolation of the other for N < 2 where the 
former breaks down. 

The first observation is that, in both cases, the max- 
imum probability fi — fc of seeing one fluxon or anti- 
fluxon is fractionally less than 50% (48.6%), occurring 
a.t N = Nc ^ 10.5. The assumption made in our JTJ 
papers is that /i (which we have called /i in our papers) 
scales linearly with N cx C/^, which we see is valid at 
best only until /i ~ 0.3. We have already just about 
achieved this in the existing experiments but have not 
yet been able to quench fast enough to provide a direct 
test of the model predictions of Fig. 4. There is, however, 
a problem that stops us embracing the Gaussian corre- 
lation approach wholeheartedly, as we have presented it 
here. As we have noted, the assumption of Gaussian 
winding number density can only be approximate, since 
it does not take periodicity (mod 27r) into account. We 
have seen that this does not matter in the calculation of 
fo and /i. 

However, for small annuli there is a problem with (|16p 
in that Fourier components f±2iN) are slightly negative 
(l/ifcal < 0.01) for very small N (with a similar problem 
for the very much smaller f^). On the other hand, by 



IV. FLUXON PRODUCTION IN AN 
EXTERNAL FIELD S / 

Let us now apply a perpendicular uniform magnetic 
field B to the A JTJ. This breaks the (j) ^ —cj) symmetry 
that is equally, the refiection symmetry of the system 
in the plane of the barrier. Once superconducting, the 
AJTJ expels the magnetic field, but we assume that a 
small fraction e of the applied field 'leaks' in the radial 
direction through the barrier, forming fluxons. The effect 
of this fleld is to produce an non-zero average winding 
number (n) = n{B). 

The result is a shift in phase gradient along the annulus 
(coordinate x) of the form: 

d^^^d^4>+'^A^, (19) 
he 

where = — A'^ is the jump in the vector potential 
across the oxide layer. 

If C = 27ri? is the circumference of the ring, radius R, 
then the change in (f> due to B is: 

Acjj^e'^ <fdl.A = ^^C^B. (20) 
he J 2TThe ^ ' 

The change in fluxon number is: 

niB)^l-A^^^^eC^B. (21) 

(For future purposes, let us call Bi the fleld value for 
which A(j> = 27r, that is h{Bi) = l.)When h is small. 
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e adjusts in any individual experiment so as to make n 
integer. As a first approximation we take the fraction e 
to be independent of B. 

In the presence of external fields we are not primarily 
interested in the small N linear regime and it is suffi- 
cient to work with Gaussian probabilities. The natural 
extension of G'($) of ^ for an AJTJ in a perpendicular 
magnetic field B is that the phase distribution will still 
be normal with variance a'^{N,n{B)) cx N, where we re- 
tain the definition of TV oc C/f of the previous section, 
but with non-zero average = 2TTfi{B): 



exp 



2a^{N, n) 



The trapping probabilities in the presence of an external 
magnetic field will then be: 



fra{Njl) 



7r+2m7r 



(22) 



— 7r+2m7r 



To a first approximation we assume that a'^{N,n) is 
independent of n{B), as would follow from assuming 
Gaussian correlations, i.e. for integer fi, fn-m{N ^n) — 
fn+jn{N,n) = f±m{N,0). This allows us to repeat the 
identification a^{N)/4:TT^ = N/12. 

As before, we are primarily interested in the proba- 
bility of seeing a single fluxon or antifluxon, but now 
for fixed tq (or N), as a function of B (or n). Now, 
as far as n 7^ 0, then the symmetry is broken and 
f+i{N,n) 7^ f^i{N,n). More precisely, in Fig. 5 we 
show /i {N, n) as a function of n for fixed N for several 
N, as derived from Gaussian distributions Eq. 
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FIG. 5: The values of fi{N, n) as a function of n for fixed A'^ 
for several values oi N; N = 1, 2, 4, 6, 8, 10, 16, according to 

Eq. (ng. 

The main characteristics of Fig. 5 are that: 



1. for N < Nc ~ 10.5 there is a double peak, cor- 
responding to the ensemble average production of 
a single fluxon by the applied external magnetic 
field. This is understood as follows. Essentially, 
when, for B = Bi we have n = 1, the zero-field 
no-trapping frequency fo{N,0) becomes the single 
fluxon trapping f+i{N,l) = /o(7V, 0), giving the 
righthand peak for positive n « 1; reversing the 
field, /_i(7V, -1) = /o(iV,0) and we get the_peak 
for negative fi « — l.The minimum value of /i be- 
tween the peaks is fi{N), as given by the KZ sce- 
nario. 

2. as N increases to Nc the height /i(iV, O)drops. 
The variance also increases, the $ distribution gets 
broader and the two peaks in n of /i {N, n) merge 
at N = Nc at the value fc ~ 0.5; 

3. as N increases beyond Nc, there is only a single 
peak centered on n = 0. We now see that the rea- 
son why the probability of seeing a fluxon decreases 
as \B\ increases is a consequence of an increased 
ability to create more than one fluxon. 

The curves in Fig. 5 for = 1 and A = 16 bear a strong 
resemblance to the experimental data shown in Fig. 2 and 
Fig. 3, respectively. They comply with the most simple 
qualitative test of our analysis, that the double peaks in 
Fig.l occur at a /ii(7/ier frequency than fc ~ 50%, and the 
single peak in Fig. 2 at a lower frequency than 50%, as 
predicted. Further, we also found that the defect winding 
number flips when we move from the left to the right 
peak, as foreseen by both models. This discrimination 
between a trapped fluxon or antifluxon can be achieved 
by measuring the transverse magnetic field dependence of 
the junction critical current Ic{B) (the details of this new 
effect and its theoretical interpretation will be reported 
elsewhere [TtI). 

More specifically, we note that the 2.0mm long sample, 
being 4 times longer than the 0.5mm long sample, ac- 
cording to (|2ip . should be 16 times more sensitive to 
the externally applied magnetic field i?, if €{B) is inde- 
pendent of B and identical for both samples. There is, 
indeed, a strong difference in sensitivity, but only by a 
factor of 7 — 8, showing that these assumptions are ap- 
proximate. What is more difficult to understand quan- 
titatively is the 16-fold increase in A^ for a four-fold in- 
crease in perimeter. This requires the efficiency factor 
a relating A^ to C to vary by a factor of 4 between the 
samples or, more fundamentally, that A^ is not linear in 
C. Of itself, the latter does not change the scaling be- 
haviour of ([2]), but the scaling exponent a. In [13, [HI 
we showed that the observed value for a was not that for 
idealised JTJs [l^. We explained this as a consequence 
of fabrication methods, but this reopens the issue. 

To go further, and match the data profiles better, 
we need specific properties of JTJs, beyond the gener- 
ics of the KZ picture (or Gaussian correlations). In 
particular, the assumption of a'^{N,h{B)) being inde- 
pendent of n{B) is oversimple. More realistically [3l| . 
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o-2(Ar, n{B)) = a'^{N, 0)(1 + fc^B^), where fc^ ex (Aj/C)^, 
i.e. it is inversely proportional to the ring area and to 
the Josephson current density Jc, since Aj cx l/Jc- Both 
are consistent with the experimental data, such as the 
flattening of /i for small B field values and the fact that 
the peak amplitudes /±„(±i?„) strongly decrease with n, 
as we shall see elsewhere. 



V. CONCLUSIONS 

We have developed two complimentary theoretical 
approaches to understand the experimentally observed 
spontaneous production of fluxons on quenching annular 
JTJs in the presence of an externally applied transverse 
symmetry breaking magnetic field B. They either as- 
sume Gaussian probabilities (as motivated by KZ causal 
horizons) or Gaussian correlation functions (as motivated 
by models for transitions based on the rapid growth of 
instabilities). Both of these approaches, which are, ap- 
proximately, identical, lead to the same scaling behaviour 
of ll]), from which 12]) follows in the appropriate regime. 

The theory is able to nicely reproduce the double peak 



behavior of the likelihood /i to produce a single de- 
fect (fluxon) shown in Fig. 2 for a sample having a circum- 
ference equal to 0.5mm and the single peak in Fig. 3 for a 
sample of circumference 2.0mm. When we began exper- 
iments on fluxon production in an external symmetry- 
breaking field, we anticipated the behaviour shown in 
Fig. 2, and not that of Fig. 3, which was initially incom- 
prehensible. We now understand it, as a consequence 
of the ease of producing more than one fluxon in larger 
annuli. Specifically, our models do provide a good first 
approximation at a better than qualitative level and the 
experimental success of the Gaussian picture in describ- 
ing the production of fluxons is of a piece with the scaling 
behaviour so robustly demonstrated in Fig.l. 
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